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Abstract 

Electronic states in silicon quantum dots are examined theoretically, taking into account a mul- 
tivalley structure of the conduction band. We find that (i) exchange interaction hardly works 
between electrons in different valleys. In consequence electrons occupy the lowest level in different 
valleys in the absence of Hund's coupling when the dot size is less than 10 nm. High-spin states 
are easily realized by applying a small magnetic field, (ii) When the dot size is much larger, the 
electron-electron interaction becomes relevant in determining the electronic states. Electrons are 
accommodated in a valley, making the highest spin, to gain the exchange energy, (iii) In the pres- 
ence of intervalley scattering, degenerate levels in different valleys are split. This could result in 
low-spin states. These spin states in multivalley artificial atoms can be observed by looking at the 
magnetic-field dependence of peak positions in the Coulomb oscillation. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 



Quantum dots fabricated on semiconductor heterostructures, e.g., GaAs/AlGaAs, are 
called artificial atoms. The electronic states in such quantum dots have been elucidated by 
examining the Coulomb oscillation; a peak structure of the current as a function of gate 
voltage which changes the number of electrons in the dots, A^, one by one. In disk-shaped 
quantum dots, discrete energy levels show an atomic-like shell structure.^ They are filled 
consecutively with increasing A^. When a shell is partly filled, electrons occupy degenerate 
energy levels with parallel spins (Hund's rule), as seen in the periodic table of atoms.^^ 
These properties of the artificial atoms have been quantitatively explained by theoretical 
studies including the electron-electron Coulomb interaction exactly .^i^i^i^i'^'i^i^ 

Compared with the electronic states in GaAs quantum dots, those in silicon (Si) quantum 
dots have been poorly understood. A combination of microfabrication and oxidation tech- 
niques on Si enables to fabricate so small quantum dots as ~ 10 nm .-^^i-^^i-^^i-^^i-^'^i^^i^^ The Si 
quantum dots are being extensively studied for their application to single-electron devices 
working at room temperature. However, the lack of information on the precise shape of 
individual dots has prevented fundamental researches of the electronic states. Particularly, 
there have been few theoretical studies which consider the effects of a multivalley structure 
of the conduction band on the electronic states in Si quantum dots.^^ The multivalley struc- 
ture should lead to different properties of electronic states from those in GaAs quantum dots 
with single valley at F point. 

Recently some experiments of the Coulomb oscillation have been reported for Si quantum 
^Q^g,i4ji5ji6 rpj^g magnetic-field dependence of the current peaks indicates various spin states 
in the quantum dots. Another experiment has implied the observation of the Kondo effect, 
which might be due to degenerate energy levels by the multivalley structure of the conduction 
band.-ifi- The explanation for these phenomena requires the understanding of the electronic 
states in Si quantum dots. The study of spin states is also important for the application of 
the quantum dots to "spintronics," quantum computers,— etc. 

The purpose of the present paper is to elucidate fundamental characters of the electronic 
states and spin configurations in Si dots, "multivalley artificial atoms," for the first time. 
We adopt the effective mass approximation,— assuming that the size of quantum dots is 
sufficiently larger than the lattice constant, to consider the following properties of Si; (i) 
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more than one valley of the conduction band (two equivalent valleys in Si-MOS systems 
and six in the bulk), (ii) anisotropy of the effective mass at a valley (m^* = 0.98mo, = 
0.19mo), and (iii) large Zeeman effect (effective g-factor is ~ 2 in Si, whereas g* ~ 0.4 in 
GaAs). We determine the electronic configuration of the ground state, considering Coulomb 
and exchange interactions. We show that the exchange interaction hardly works between 
electrons in different valleys. In consequence electrons occupy the lowest level in different 
valleys in the absence of Hund's coupling when the dot size is less than 10 nm. High-spin 
states are easily realized by applying a small magnetic field. This is quite different from 
the case of single-valley artificial atoms of GaAs quantum dots, where high spins appear 
due to the exchange interaction only in the case of orbital degeneracy (Hund's rule). When 
the dot size is much larger than 10 nm, the electron-electron interaction becomes relevant 
in determining the electronic states. Electrons are accommodated in a valley, making the 
highest spin, to gain the exchange energy. 

We also investigate the effects of anisotropic shape of quantum dots and intervalley 
scattering. The former tends to make higher spins, whereas the latter could make lower 
spins. These spin states in multivalley artificial atoms should be observed by looking at the 
magnetic-field dependence of the peak positions in the Coulomb oscillation. We calculate the 
addition energies, as functions of the magnetic field, which correspond to the peak positions 
in the Coulomb oscillation. 

The organization of the present paper is as follows. In the next section (Sec. |H)), we 
explain our model for Si quantum dots and calculation method using the effective mass 
approximation. In Sec. IIIH we present the calculated results for the electronic states in 
isotropic Si quantum dots. We discuss the fundamental properties of multivalley artificial 
atoms, with changing the size of the dots. In Sec. IIVI we examine anisotropic Si quantum 
dots. In Sec. El effects of intervalley scattering on the electronic states are considered. The 
last section (Sec. IVIjl is devoted to our conclusions and discussion. 

II. MODEL AND CALCULATION METHOD 

For the multivalley structure of conduction band in Si, we assume that two valleys at k = 
(0, 0, ±/co), where ko = 0.85 x 27r/a (a is the lattice constant, 0.543 nm), are equivalent and 
located below the other valleys at k = (ifco, 0, 0) and (0, ifco, 0) by Ai^vaiioy Ai^vaiicy = 100 
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meV, considering the effect of oxidation-induced strainjii*2i The valleys of = (0,0,±fco), 
(i/co, 0, 0), and (0, ±A;o, 0) are denoted by ±kz, ±kx, and ±ky hereafter. 

As a model for the confinement potential of quantum dots, we adopt a three-dimensional 
harmonic potential V{x,y, z). We examine an isotropic confinement 



(1) 



to focus on generic properties of electronic states in Si dots, for a while. An anisotropic 
confinement is studied later (Sec. llVjl . 

First, we calculate one-electron states and energy levels in the absence of electron-electron 
interaction. In the effective mass approximation,— the wavefunction for electrons around a 
valley, e.g., ±kz, is described by a product of the Bloch wave at the valley and the envelope 
function F-j.^^ 



(2) 



where = (0,0, /cq) and u±k^{r) is a function with the periodicity of the Bravais lattice. 
The envelope function F±kz satisfies the effective mass equation 

+ 



2m* dx"^ 



dy'^ 2m^ dz"' 



+ V{x,y,z) 



(3) 



where the energy e is measured from the bottom of conduction band at ±A;^. [In the 
effective mass equation around the valley of -^k^ {±ky), rrii appears in the differential term 
with respect to x {y)-] We disregard the interband scattering on the assumption that the 
confinement potential V{x,y, z) is so smooth that its Fourier components with A; ~ 1/a are 
negligibly small. (The effect of interband scattering, which may stem from impurities within 
the dots, sharp edges of the dots, etc., is discussed in Sec. |Vl) Using V{x,y,z) in Eq. ((H), 
Eq. Q yields the eigenvalues of 
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with Tlx, riy, Uz — 0, 1, 2, • • • . [For valleys of ±kx {^ky), m* appears in the first (second) 
term.] The eigenfunction for the lowest level with {ux, riy, Uz) — (0, 0, 0) is given by 



X exp 



(5) 

The size of the confinement potential is given by 



I = 2^Jh/mluu (6) 

where ujt — yjKjmX^ using the smaller effective mass m^. Three cases of I =5 nm, 10 nm, 
and 15 nm are examined. 

Many-body states are simply represented by single configuration of the occupation of 

the one-electron levels by electrons, for semi-quantitative discussion. The Coulomb and 
exchange interactions between electrons {ki;nx,ny,nz) and {kj;n'^,ny,n'^) are written as 





(7) 



and 



1 1 ^^ndr2C"""^*(n)V'r;'""""(^2) 
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Ane I ri — r2 I 

X F;;'"--(r2) uUr,)ul^{r,)uu,{r,)u,^{r,\ (8) 

respectively, with dielectric constant e — II.Osq- These terms are evaluated numerically. We 
calculate the total energy for all possible configurations and determine the ground state. 
In a magnetic field B = (0, 0, B), we consider the Zeeman energy 

-Ezeeman = —g*fJ'BStot,zB, (9) 



where (?* = 2, /iB is the Bohr magneton, and Stot,z is the z component of the total spin. The 
orbital magnetization effect is disregarded because of large effective mass in Si.— 

III. ELECTRONIC STATES IN SI QUANTUM DOTS 

In this section, we clarify the basic properties of electronic states in Si quantum dots, 
considering an isotropic confinement, Eq. (|T|). 

A. Electronic states 

We begin with the one-electron energy levels. As shown in Eq. (0)), the one-electron levels 
are labeled by valley index and quantum numbers (n^., n^^, n^) for the orbital motion. The 
lowest levels belong to the valleys of +kz or —k^ with {nx,ny,nz) = (0,0,0) (see horizontal 
lines in Fig. ^). The next levels are (itA;^; 0, 0, 1), which are lower than levels (±/c^; 1, 0, 0) 
and (ifc^; 0, 1, 0) in spite of the isotropic confinement. This is due to the anisotropy of the 
effective mass [e{±k,; 1, 0, 0)-e{±k,; 0, 0, 1) = h{^K/m*- ^K/m*)]. We find that all the 
electrons are accommodated in valleys of ±A;^ for = 1 to 5 in all the cases of the dot size. 

Regarding the electron-electron interaction between different valleys, we obtain two sig- 
nificant results, (i) The Coulomb interaction between {±kz',nx,ny,nz) and {±kz;n'^,ny,n'^) 
is the same as that between {±kz]nx,ny,nz) and {^kz;n'^,n'y,n'^). This is because the en- 
velope functions are equivalent for the valleys of +kz and —k^ in Eq. (|7j). (ii) The exchange 
interaction is negligibly small between electrons in different valleys, owing to the integra- 
tion of rapidly oscillating factors in Eq. (jSI). Therefore, there is no spin coupling (Hund's 
coupling) between different valleys. 

When the dot size is Z = 5 nm or 10 nm, two electrons occupy the lowest level in different 
valleys in the absence of Hund's coupling. In consequence, the energies of different spin 
states are degenerate. The configuration in Fig. [Tfai) has the total spin of Stot = 0, whereas 
that in Fig. ^a2) has Stot = l/2cg)l/2 = 0©l. These configurations have the same energy. 
In a small magnetic field, the Zeeman effect, Eq. 0, makes the largest spin of Stot = 1 with 

>S'tot,z = 1- 

When the dot size is / = 15 nm, the ground state with N = 2 has another configuration; 
one electron occupies an upper level to gain the exchange energy with the other electron 
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with parallel spin in a valley [Fig. ^b)]. Note that the exchange interaction, Eq. (jSJ, works 
effectively between electrons only in the same valley. In general, the electronic states are 
determined by the competition between the electron-electron interaction and spacing of one- 
electron energy levels. When the dot size is /, the former is characterized by E^^t = e^/ {Airel), 
whereas the latter is by -Eieveis = h'^/{m^l^). The ratio of these energies, -E'int/-E'ioveis = 
mj'e^//(47re^^), is 1.52, 3.04, and 4.56 for / = 5 nm, 10 nm, and 15 nm, respectively. When 
the ratio is larger, more electrons are accommodated in a valley with parallel spins to gain 
the exchange energy. 

The ground state with = 3, shown in Fig.|2l can be understood in a similar way. For the 
smallest dot (/ = 5 nm), three electrons occupy the lowest levels [S'tot = 1/2; Fig. |2fa)]. For 
the largest dot (/ = 15 nm), they occupy upper levels in a valley, making S'tot = 3/2, to gain 
the exchange energy [Fig.|21^b)]. For the intermediate case (/ = 10 nm), two electrons make a 
spin triplet, occupying two levels in a valley, while an electron is accommodated in the lowest 
level in the same valley (ci) or in the other equivalent valley (C2). These configurations have 
the same energy but different spin, (ci) S'tot = 1/2 and (C2) S'tot = l(S>l/2 = l/2©3/2, the 
largest spin of which is realized in a magnetic field. 

With increasing magnetic field, the Zeeman effect, Eq. (jH)), makes larger spin states 
more favorable. For = 3, the ground state could change from configuration Fig. Efa) 
(•S'tot = 1/2) at B < Be to configuration Fig. |2fb) (S'tot = 3/2) at i? > Be- Such a transition 
is observed at Be = 9.47 T when the dot size is / = 7.5 nm. (The same transition takes 
place at much larger Be when / = 5 nm.) We find a similar transition with = 4, from 
S'tot = 1 to S'tot = 2 at i?c = 3.44 T, when the dot size is / = 7.5 nm. 

B. Addition energies 

For the comparison with the experimental results, we calculate addition energies which 
correspond to peak positions of the Coulomb oscillation.ii^ The addition energy of the iVth 
electron on the dot is given by 

Hn = En — E^-i, (10) 
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where En is the ground-state energy with electrons. Figure El shows fiN (A^ = 1 to 5), as 
functions of magnetic field B. The slope of iin{B) is given by 

-g*^iB[Stot{N)-S,^,{N -1)1 

which reflects the total spin S'tot(^) with N electrons. 

In Fig. |2fa), the dot size is the smallest (/ = 5 nm). The one-electron levels are consec- 
utively occupied with increasing the number of electrons A^. For N = 2, the electrons are 
accommodated in different valleys with parallel spins in a magnetic field, as discussed previ- 
ously. In consequence, Stot{N) = 1/2, 1, 1/2, 0, and 1/2 for = 1 to 5, respectively. This 
results in the slope of —g*fiB/'2 for /ii(-B), /i2(-B), fJ^^^B), and +g*fiB/2 for fis^B), fi4{B). 

In Fig.El^b), the dot is the largest (/ = 15 nm). To gain the exchange energy, the electrons 
occupy upper levels in a valley, making the highest spin, Stot{N) = N/2. Accordingly, iin{B) 
[N = 1 to 5) are parallel to each other, with the same slope of —g*fiB/'2- When the dot 
size is / = 10 nm, the highest spin of Stot{N) = N/2 is also realized. The slopes of fJ,N{B) 
{N = 1 to 5) are the same as those in Fig. Efb). 

In Fig. Ofc), the dot size is / = 7.5 nm. fisiB) shows a kink at B^ = 9.47 T, which is 
due to the transition of the ground state with A^ = 3 from 5'tot(3) = 1/2 to 3/2. Since 
5^1(2) = 1, the slope of fJ^s^B) is g*fiB/'^ and —g*fiB/2 at B < B^ and B > B^, respectively. 
This transition results in a kink of ^a{B) at the same magnetic field [with S'tot(4) = 2, the 
slope of ^ii{B) is — 3(7*/iB/2 at i? < B^, —g*fiB/'2 at B > Be]. Similarly, the transition of the 
ground state with A^ = 4 makes a pair of kinks on fi4{B) and /i5(-B) at B^ = 3.44 T. Such 
pairs of kinks of adjacent addition energies have been observed in the Coulomb oscillation.— 

At 3.44 T< B < 9.47 T in Fig. Efc), fi^iB) has a slope of -3g*HB/2 (indicated by 
dotted line). Then the conductance through the quantum dot should be suppressed since 
the addition of the fourth electron is forbidden by the spin selection rule: S'tot = 1/2 with 
A^ = 3 and S'tot = 2 with A^ = 4. This is called spin blockade.-^ The spin blockade in Si 
quantum dots has been found experimentally.— 

IV. ANISOTROPIC CONFINEMNET 

In the previous section, we have examined the isotropic confinement to illustrate the 
basic properties of Si quantum dots. In this section, we investigate an anisotropic con- 
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finement. Quantum dots, which are fabricated in the middle of Si wires with oxidation 
technique /°i^^i^^i^^i^^i^^ may be elongated along the wires. We adopt an elliptical potential, 
elongated in [110] direction of the lattice,— 



V{x,y,z) 



K' 



X + y 



+ K 



-X + y 



(11) 



with K' < K. By the transformation of the axes to X = (x + y)/V2, Y = ( —X + y)/ v2, 
and z, it is rewritten as V{X,Y,z) = {1/2){K'X^ + KV^ + Kz"^). For valleys of ±k^, the 
one-electron energy levels are given by 



e{±k^;nx,nY,nz 




K' 1 
h^|—{nx + -) + h^ 
mt 2 



mJ 2 




^ ,K . 1, 



(12) 



We change the size of the confinement potential in [110] direction, Ix = 2^ hjmXijJxt where 



LiJxt = ^/K' /m^, as 5 nm, 10 nm, and 15 nm. The size in [110] and [001] directions is fixed at 
ly = Iz = ^ nm, where ly = 2^Jh/m*ujYt) and Iz = 2\Jh/mluJzt with uyt = ^zt = a/ K/m^. 

Figure 0] shows the electronic configuration of the ground state with N = 4 when (a) 
Ix = 5 nm (isotropic dot), (b) Ix = 10 nm, and (c) 15 nm. As the shape of the dot is 
more elongated in [110] direction, the energy levels (nx,0,0) {nx = 1,2,3. . .) are lowered 
[see Eq. (fT^ with smaller K']. Electrons begin to occupy these levels: (a) (iA;^; 0,0,0) 
only, (b) (iA;^; 0,0,0) and (i/^z; 1, 0, 0), and (c) (0,0,0) to (3,0,0) in a valley. This means 
that the electrons tend to be populated in the elongated direction. Reflecting small level 
spacings between (nx,0,0), high-spin states are realized more often: (a) S'tot = 0, (b) 
Stot = 1C>?)1 = 0©1©2, and (c) S'tot = 2. We find that the ground state has the largest spin 
of Stot{N) = N/2 for X = 1 to 5 in elongated dots with Ix = 10 nm and 15 nm [ly = Iz = 5 
nm) in a magnetic field. The magnetic-field dependence of addition energies shows the 
same slope of —g*^^/2, which is similar to the situation in Fig. Efb). 



V. INTERVALLEY SCATTERING 



In the present section, we discuss the effects of intervalley scattering. The intervalley 
scattering is not relevant when the size of the dots is much larger than the lattice constant 
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and confinement potential is sufficiently smooth. However, the scattering by impurities 
within the dots, sharp edges of the dots, etc., mixes the electronic states in different valleys. 
As an example, let us consider an impurity potential at R inside a quantum dot, 

V,^p{r) = Vo6{r-R), 

to the lowest order of the perturbation. The degenerate levels of {±kz',nx,ny,nz) are split 
into two by this potential, 

£:(a; nx, Uy, n^) = ^(ibfc^; Ux, Uy, n^), 

and 

s(b, Tlx, ^z) ^(i^2) '^xi '^yi '^z) ~l~ ^nx,ny,nz i 

where A^^^^^^^^ = 2Vo\F^j^^^'^^ {R)\'^\u±kXR)\'^ ^ as shown in Fig. Efa). The corresponding 
wavefunctions are expressed as 

with 6 = a.Tg{+kz; Ux, Uy, n^ll^mpl — kz] Ux, Uy, n^). For the electron-electron interactions, we 
find that (i) the Coulomb interaction between {a;nx,ny,nz) and {a;n'x,ny,n'^) is identical 
to that between (6; Ux, riy, n^) and (6; n'^, n'y, n'^), and also to that between (a; Ux, riy, Uz) and 
{b;n'^,ny,n'J. (ii) The exchange interaction is negligibly small between {a;nx,ny,nz) and 
{b;n'x,n'y,n'^). 

The intervalley scattering influences the spin states in a spherical dot of Z = 5 nm. In 
the absence of intervalley scattering, the first two electrons occupy the degenerate levels 
of (±/c^; 0,0,0) with parallel spins [Stoti'^) = 1] in a small magnetic field [Fig. ma2)]. In 
the presence of intervalley scattering, the two electrons occupy the lowest level of (a; 0, 0, 0) 
with anti-parallel spins [Stot{2) = 0], as shown in Fig. El^a). For = 1 to 5, the spin states 
become S'tot(A^) = 1/2, 0, 1/2, 0, and 1/2, respectively. 

The addition energies are shown in Fig. Efb), as functions of magnetic field, in the case 
of Ao,o,o = 2 meV. The slopes of fiN are —g*fiB/'2 or +g*fiB/2 alternately with increasing N, 
reflecting the above-mentioned spin states. In the case of Ao,o,o = 0.4 meV, the transition 
of the ground state with = 2 is observed at the magnetic field of = 3.45 T, where 
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g*jji^Bc = Ao,o,o- For B < Be, two electrons occupy level (a; 0,0,0), making 5*101(2) = 0, as 
shown in Fig. |Sl^a). For B > Be, one electron occupies level (a; 0, 0, 0) and the other occupies 
level (6; 0,0,0) with parallel spins [S'tot(2) = 1], to gain the Zeeman energy. This transition 
makes a pair of kinks on /i2 and fi^ [see Fig. EI^c)]. 

When the dot size is much larger than 10 nm, the intervalley scattering never changes the 
spin states. Electrons occupy the levels of (a; 0, 0, 0), (a; 0, 0, 1), (a; 1, 0, 0), ■ ■ ■ , consecutively 
with parallel spins to gain the exchange energy [not the levels of (6; n^, Uy, Uz) owing to the 
absence of exchange interactions between (6; rix, Uy, n^) and (a; n'^, n'y, n'J]. Then the largest 
spins of Stot{N) = N/2 are realized. The addition energies Hn^B) (A^ = 1 to 5) have the 
same slope of — 5'*/iB/2, in the same way as in Fig. Efb). 

VI. CONCLUSIONS AND DISCUSSION 

We have theoretically examined the electronic states in Si dots, multivalley artificial 
atoms, (i) In spherical quantum dots smaller than ~ 10 nm, the one-electron levels in 
equivalent valleys are occupied consecutively with increasing the number of electrons A^. 
High-spin states are attributable to the absence of spin couplings between different valleys. 
This is quite different from the case of GaAs quantum dots, where high spins appear due to 
the exchange interaction in exceptional cases with orbital degeneracy (Hund's rule), (ii) In 
quantum dots much larger than 10 nm, the electrons tend to occupy upper levels in a valley, 
making the largest spin of ^tot = N/2, to gain the exchange energy, (iii) In elliptical quantum 
dots, the electrons tend to be populated in the elongated direction. High-spin states often 
appear, (iv) In the presence of intervalley scattering, the degenerate energy levels in different 
valleys are split. In small dots, this results in the lowest spin states (^tot = or 1/2). The 
spin states in large dots, S'tot = N/2, are not influenced by the intervalley scattering. 

In the present study, the electronic states are represented by single configurations of 
the occupation of one-electron levels by electrons, for semi-quantitative discussion. We 
have neglected the configuration interaction, or correlation effect. (It should be noted that 
the configuration interaction vanishes between electrons in different valleys for the same 
reason as the exchange interaction vanishes.) In this approximation, the energy of low-spin 
states is overestimated compared with the energy of high-spin states. Hence we should have 
underestimated the dot size where the highest spins appear. The correlation effect becomes 



11 



more important as the dot size is larger (E'int/i^ieveis is larger) for quantitative discussion.— 
The consideration of the correlation effect with multivalley structure in Si dots requires 
further study. 

The spin states in Si quantum dots depend on both their size and shape. Our cal- 
culated results cannot be compared directly with experimental results of the Coulomb 
oscillatio n^^i^^i^^ since the exact shape of quantum dots is not known. Roughly speaking, the 
magnetic-field dependence of peak positions in the observed Coulomb oscillation is similar 
to that in Figs. |Sfb), (c) rather than that in Figs. ini^a)-(c), which might imply that low-spin 
states are realized in small dots and in the presence of intervalley scattering. Qualitatively, 
pairs of kinks have been observed,— which indicate the transitions of the ground state. The 
spin blockade with l^totl^^ + 1) " 'S'tot(^)| > 3/2 has also been reported.JA We have shown 
that the spin states can be changed by the intervalley scattering. If the strength of the 
intervalley scattering were tuned, e.g., using STM tips, the spin states in Si dots could be 
controlled. 

Recently, quantum computing devices have been proposed, utilizing electron spins in 
quantum dots.^^ For the devices, the decoherence time of electron spins must be much longer 
than the gate operation time. Main origins of the decoherence are spin-orbit interaction and 
hyperfine interaction with nuclear spins in the dots. In Si, the spin-orbit interaction is much 
weaker than in GaAs. The hyperfine interaction is also weaker due to smaller isotope ratio 
with nuclear spins (^^Si of 4.7% only). Hence a longer decoherence time is expected in Si 
quantum dots than in GaAs quantum dots. Therefore, Si dots may be advantageous to the 
quantum computing devices. In Si dots, high-spin states are realized easily, which could be 
utilized for the devices. 
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FIG. 1: Electronic configurations for the ground state with N = 2 in isotropic quantum dots. 
{nx,ny,nz) are the orbital quantum numbers for each one-electron level, whereas ±kz are the 
valley indices. Configuration (ai) or (a2) is realized when the dot size is / = 5 and 10 nm. 
They are degenerate in energy in the absence of magnetic field. The total spin is ^tot = and 
'S'tot = 1/2 <8) 1/2 = © 1, respectively. Configuration (b) with Stot = 1 is realized for I = 15 nm. 
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FIG. 2: Electronic configurations for the ground state with = 3 in isotropic quantum dots. 
{nx,ny,nz) are the orbital quantum numbers for each one-electron level, whereas ib/c^ are the 
valley indices. Configurations (a) and (b) are realized when the dot size is / = 5 nm and 15 nm, 
respectively. Configuration (ci) or (02) appears for I = 10 nm. When the dot size is / =7.5 nm, 
configuration (a) appears at i? < 9.47 T, whereas (C2) appears B > 9.47 T. 
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FIG. 3: Addition energies, {N = 1 to 5), as functions of magnetic field B, in isotropic quantum 
dots. The dot size is (a) / = 5 nm, (b) 15 nm, and (c) 7.5 nm. In (c), the spin blockade occurs for 
fi4{B) at 3.44 T<B < 9.47 T (dotted line). 
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FIG. 4: Electronic configurations for the ground state with iV = 4 in eUiptical quantum dots. 
{nx,nY,nz) are the orbital quantum numbers for each one-electron level, whereas ±kz are the 
valley indices. The size of the dot in [110] direction is (a) Ix = ^ nm, (b) 10 nm, and (c) 15 nm, 
whereas the size in the other directions is fixed at ly = Iz = ^ nm. 
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FIG. 5: (a) The effect of intervalley scattering on the one-electron energy levels. The degenerate 
levels of {±kz;nx,ny,nz) are split into two. (b), (c) Addition energies, fiN = 1 to 5), as 
functions of magnetic field B, in isotropic quantum dots of / = 5 nm. The energy splitting by the 
intervalley scattering is (b) Ao,o,o = 2 meV and (c) Ao,o,o = 0.4 meV. 
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